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Abstract 

We present an example of a scenario of particle production and decay in supersymmetry 
models in which the supersymmetry breaking is transmitted to the observable world via grav- 
itational interactions. The case is chosen so that there is a large production of tau leptons in 
the final state. It is characteristic of large tan/? in that decays into muons and electrons may 
be suppressed. It is shown that hadronic tau decays can be used to reconstruct final states. 
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1 Introduction 



If supersymmetry (SUSY) exists at the electroweak scale, then gluinos and squarks will be copiously 
produced in pairs at the LHC and will decay via cascades involving other SUSY particles. The 
characteristics of these decays and hence of the signals that will be observed and the measurements 
that will be made depend upon the actual SUSY model and in particular on the pattern of super- 
symmetry breaking. Previous, detailed studies of signals for SUSY at the LHC ||, |J] have 
used the SUGRA model ||, |(| , in which the supersymmetry breaking is transmitted to the sector 
of the theory containing the Standard Model particles and their superpartners via gravitational 
interactions. The minimal version of this model has just four parameters plus a sign. The lightest 
supersymmetric particle (xi) has a mass of order 100 GeV, is stable, is produced in the decay of 
every other supersymmetric particle and is neutral and therefore escapes the detector. The strong 
production cross sections and the characteristic signals of events with multiple jets plus missing 
energy Jfix or with like-sign dileptons £ ± ^ ± plus $t 01 enable SUSY to be extracted trivially from 
Standard Model backgrounds. Characteristic signals were identified that can be exploited to de- 
termine, with great precision, the fundamental parameters of these model over the whole of its 
parameter space. Variants of this model where R-Parity is broken [|| and the Xi is unstable have 
also been discussed §. 

These models have characteristic final states depending upon their parameters. The next to 
lightest neutral gaugino x® is produced in the decays of squarks and gluinos which themselves may 
be produced copiously at the LHC. The decay of X2 then provides a tag from which the detailed 
analysis of supersymmetric events can begin. The dominant decay is usually either ~^ hxi or 
X2 — ► i + i~X\-> which can proceed directly or via the two step decay X2 ~~ > — ► l + l~Xi- The 
latter leads to events with isolated leptons. Both of these characteristic features have been explored 
in some detail in previous studies || ||], Q . 

In the previous cases the smuon, selectron and stau were essentially degenerate. At larger values 
of tan P, this degeneracy is lifted and the f± becomes the lightest slepton. If m 1 / 2 is small enough, 
then the two-body decays X2 ~^ Xih, Xi% wm n °t be allowed, and if mo is large enough, then 
X2 ~^ ^R^- wm a ^ so n °t be allowed. Then for large enough tan/3 the only allowed two-body decays 
are X2 ~^ r± '? =F — ► t +t ~Xi- I n such cases, tau decays are dominant, and final states involving tau's 
must be used. 

The simulation in this paper is based on the implementation of the minimal SUGRA model in 



ISAJET |g . We use m = m 1/2 = 200 GeV, tan/3 = 45, A = and sgn/i = -1. The mass 
spectrum for this case is shown in Table |]. The only allowed two-body decay of X2 ^ s ^i r > so 
it has a branching ratio of more than 99%. 

The total production cross-section for this model is 99 pb at the LHC. The rates are dominated 
by the production of gg and gq final states. Interesting decays include: 



• BR{x° 2 -> rfy) = 99.9%; BR(x{ -> v T f x ) = 99. i 
. BR(xs -» X° 2 Z) =13%; BR( X3 -» tt x ) = 21%; 

• BR(g -» bbx) = 55%; BR{g -» bb 2 ) = 10%; 

• BR(g -» q R qi) = 3%; BR(g q r q r ) = 5.7%; 

• BR{q~ L - xlq) = 30%; BR(q~ R j$q) = 97%; 
Here q refers to a light quark. 
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Table 1: Masses of the superpartners, in GeV, for the case being studied. Note that the first and 
second generation squarks and sleptons are degenerate and so are not listed separately. 
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All the analyses presented here are based on ISAJET 7.37 [ 10 1 and a simple detector simulation. 
600K signal events were generated which would correspond to 6 fb _1 of integrated luminosity. 
The Standard Model background samples contained 250K events for each of tt, WZ with W — > 
ev, fiu, tv, and Zj with Z — > vu, tt, and 5000K QCD jets (including g, u, d, s, c, and b) divided 
among five bins covering 50 < Pj- < 2400 GeV. Fluctuations on the histograms reflect the generated 
statistics. On many of the plots that we show, very few Standard Model background events survive 
the cuts and the corresponding fluctuations are large, but in all cases we can be confident that 
the signal is much larger than the residual background. The cuts that we chosen have not been 
optimized, but rather have been chosen to get background free samples. 



The detector response is parameterized by Gaussian resolutions characteristic of the ATLAS [11] 
detector without any tails. All energy and momenta are measured in GeV. In the central region 
of rapidity we take separate resolutions for the electromagnetic (EMCAL) and hadronic (HCAL) 
calorimeters, while the forward region uses a common resolution: 

EMCAL W%/VE®1%, \r}\ <3 

HCAL 5O%/\/£0 3%, \r)\ <3 

FCAL 100%/v A Ee7%, \r)\ >3 

A uniform segmentation Arj = A(fi = 0.1 is used with no transverse shower spreading. Both 
ATLAS Q and CMS @ have finer segmentation over most of the rapidity range, but the neglect 



of shower spreading is unrealistic, especially for the forward calorimeter. Missing transverse energy 
is calculated by taking the magnitude of the vector sum of the transverse energy deposited in in 
the calorimeter cells. An oversimplified parameterization of the muon momentum resolution of the 
ATLAS detector — including a both the inner tracker and the muon system measurements — is 
used, viz 



SP T /P T = .y/0.016 2 + (O.OOllPr) 2 
For electrons we use a momentum resolution obtained by combining the electromagnetic calorimeter 
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resolution above with a tracking resolution of the form 



5P T /P T = (3 ZT^a ) vVo004P T ) 2 + 0.0001 
This provides a slight improvement over the calorimeter alone. 



Jets are found using GET JET [1C] with a simple fixed-cone algorithm. The jet multiplicity in 



SUSY events is rather large, so we will use a cone size of 



R= y (A??) 2 + (A(ft) 2 = 0.4 

unless otherwise stated. Jets are required to have at least Pt > 20GeV; more stringent cuts are 
often used. All leptons are required to be isolated and have some minimum Pt and \rj\ < 2.5, 
consistent with the coverage of the central tracker and muon system. An isolation requirement 
that no more than 10 GeV of additional Ex be present in a cone of radius R = 0.2 around the 
lepton is used to reject leptons from 6-jets and c-jets. In addition to these kinematic cuts a lepton 



(e or /j,) efficiency of 90% and a 6-tagging efficiency of 60% is assumed fll! ]. 

As t's are a crucial part of this analysis, they require special treatment. We concentrate on 
hadronic tau decays, since for leptonic decays the origin of the lepton is not clear and the visible 
lepton in general carries only a small fraction of the true tau momentum. Using the fast simulation, 
we first identify the hadronic taus by searching the reconstructed jet list for jets with Pt > 20 GeV 
and 1 7/| < 2.5. We then compare these jets with the generated tau momenta and assign them to a 
reconstructed tau if E T > 0.8Ej et and the center of the jet and the tau are separated by AR < 0.4. 



We then rely on a full simulation IS] of Z + jet events with Z —* tt. Events were generated 
with PYTHIA g4j and passed through the ATLAS GEANT simulation (DICE) and reconstruction 
(ATRECON) programs Jig] . The charged particles were reconstructed with the tracking and the 
photons with the calorimeter. Cuts were then applied to the invariant mass and isolation of the 
reconstructed taus. These cuts produce a rejection factor against QCD jets of a factor of 15 and 
accept 41% of the hadronic tau decays. We apply these results to the hadronic tau's identified 
in our fast simulation on a probabilistic basis. The accepted hadronic decays are assumed to be 
measured using the resolution from the full simulation, while the ones not accepted are put back 
into the jet list. Fake t's are made by reassigning jets with the appropriate probability. The full 
simulation also indicates that the tau charge is correctly identified 92% of the time. We include this 
factor in our fast r reconstruction and assign the fake tau's to either sign with equal probability. 
For cases where the tt invariant mass is to be measured, the generated tt invariant mass is smeared 
with a resolution derived from the full simulation, i.e., a Gaussian with a peak at M = 0.66M rr 
and a/M = 0.12. In cases where the measured momentum of the r decay products is needed, the 
measured jet energy is used. 

Results are presented for an integrated luminosity of 10fb _1 , corresponding to one year of 
running at 10 33 cm~ 2 s ; pile up has not been included. We will occasionally comment on the 
cases where the full design luminosity of the LHC, i.e. 10 34 cm _2 s _1 , will be needed to complete 
the studies. For many of the histograms shown, a single event can give rise to more than one entry 
due to different possible combinations. When this occurs, all combinations are included. 

The rest of this paper is organized as follows. We first illustrate how measurement of the tt 
final state can be used to infer information on the masses of the staus. We then use this final 
state in conjunction with 6-jets to reconstruct gluinos and bottom squarks. Methods for extracting 
information on light squarks are then shown and the dilepton mass distribution is used to give 
information on the masses xt- Finally we show how this information can be combined to constrain 
the underlying model parameters. 
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2 Effective mass distribution 

The first step in the search for new physics is to discover a deviation from the Standard Model 
and to estimate the mass scale associated with it. SUSY production at the LHC is dominated 
by gluinos and squarks, which decay into multiple jets plus missing energy. A variable which is 
sensitive to inclusive gluino and squark decays is the effective mass M e g, defined as the scalar sum 
of the -Pr's of the four hardest jets and the missing transverse energy fa*, 

M eS =pT,l+ PT,2 + PT,3 + PT£ + fa ■ 

Here the jet Pt's have been ordered such that pr,i is the transverse momentum of the leading jet. 
The Standard Model backgrounds tend to have smaller fa, fewer jets and a lower jet multiplicity. 
In addition, since a major source of fa is weak decays, large fa events in the Standard Model tend 
to have the missing energy associated with leptons. To suppress these backgrounds, the following 
cuts were made: 

• fa > 100 GeV; 

• > 4 jets with P T > 50 GeV and p T ,i > 100 GeV; 

• Transverse sphericity St > 0.2; 

• No fjL or isolated e with Pt > 20 GeV and \rj\ < 2.5; 

• fa > 0.2M cS . 

Note that some of these jets could result from hadronic tau decays. With these cuts and the 
idealized detector assumed here, the signal is much larger than the Standard Model backgrounds 
for large M e g, as is illustrated in Figure |]. Thus, the discovery strategy developed for low tan/3 j| 
also works for this case. As demonstrated in more detail elsewhere [|| the shape of this effective 
mass distribution can be used to estimate the masses of the SUSY particles that are most copiously 
produced; here quarks and gluinos. 

3 Tau- Tau invariant mass 

As can be seen from the decays listed above we expect significant prodution of X2 an d hence of tau 
pairs from the decay of q~L We require that the events contain at least two jets that are indentified 
as hadronic tau decays using the above algorithm. In addition, the following cuts are applied: 

• fa > 100 GeV; 

• at least four jets with Pt > 50 GeV and at least one jet pr,i > 100 GeV; 

• M eS > 500 GeV; 

• fa > 0.2M cS . 

Again, some of these jets could result from hadronic tau decays. 

We then search for taus that decay hadronically using the algrotihm discussed above. The 
reconstructed rr invarient mass distribution is shown in Figure g all combinations of tau charges 
are shown in this Figure. It can be seen from this distribution that there is a clear structure. 
There is considerable background from combinations where one of the identified tau jets is from 
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Figure 1: SUSY signal (open histogram) and Standard Model backgrounds (filled histogram). 
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Figure 2: Reconstructed rr mass distribution. All combinations of tau pairs are shown irrespective 
of the charge. The dashed histogram shows the comibination of one real tau and one fake tau. The 
actual experiement would observe the sum of the two histograms. The background from Standard 
Model processes is very small and is not shown. 

a tau and the other is from a misidentified jet. The invarient mass distribution of these pairs is 
also shown in Figure it is rather featureless. The tau algorthm has not been optimized so this 
background could well have been overestimated. The background from events where both taus are 
misidentified jets and the Standard Model background are both negligible. The position of the 
peak in this mass distribution enables one to infer the position of the end point arising from the 
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Figure 3: Reconstructed rr mass distribution showing the effect of rescaling the generated tau-tau 
invarient mass distribution by ±7.5%. 
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Figure 4: Reconstructed r + r~ — r ± T ± mass distribution. The dashed line shows the fake-real 
background. The flucuations are slightly larger than the true statistics. 



decay chain \2 — ► — * Xi rr: 



M™ ax = M,o 

In order to estimate the precision with which this endpoint can be determined, the generated 
tau-tau invarient mass distribution was shifted by ±7.5% from its nominal value. The effect on the 
reconstructed rr mass distribution is shown in Figure U. These cases can clearly be distinguished. 
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Figure 5: Reconstructed tt +jet mass distribution where the jet is tagged as a 6-jet. The background 
from Standard Model processes is negligible. 



The actual precision that can be obtained on the position of this end point requires a more detailed 
study. Tau decays are well understood; the problem is to determine the effects of the detector 
resolutions and the cuts. For the purposes of extracting parameters below, we will assume an 
uncertainty of 5% can be achieved. 

There are some events beyond this edge as can be seen by looking at the subtracted distribution 
t + t - — t~t~ — t + t + shown in Figure ||. This subtraction also elliminates the background from 
fake taus because their charges are not correlated. Here the excess extends to ~ 150 GeV and is 
due to X3 and Xi decays. This can be confirmed by the large Z signal (see below). The fluctuations 
in this histogram reflect the generated statistics, which correspond to about 6fb _1 ; three years at 
low luminosity would make this high-mass signal much clearer. 



4 Reconstruction of g — > bb — > bx%b — > bbr + r Xi 

The event sample of the previous section is used in an attempt to reconstruct squarks and gluinos. 
We concentrate here on final states with b quarks as these have the larger branching ratios and less 
combinatorial background. In addition to the previous cuts, we require a tagged 6-jet with Pt > 25 
GeV; this jet could be one of the ones in the previous selection. Events are selected that have 
reconstructed tau pairs with invariant mass within ±10 GeV of peak in Figure ^, and the invariant 
mass of the tau pair and the 6-jet is formed. This mass distribution is shown in Figure |. The 
sign subtracted distribution corresponding to t + t~~ —t~t~ — t + t + is used to reduce combinatorial 
background. There should should be an edge at ~ — = 310 GeV. The edge is not sharp 
— 3 particles are lost, two u T 's and the Xi- I n addition the distribution is contaminated by decays 
from X3 an d Xa- The structure is not clear, but is well distinguished from that resulting from the 
case where the 6-jet is replaced by a light quark jet, shown in Figure || 

Further information can be obtained by applying a partial reconstruction technique. This was 
developed in Ref. [l] (so called "Point 3") where the decay chain g — > bb — > bbx^b — > bb£ + £~Xi was 
fully reconstructed as follows. If the mass of the lepton pair is near its maximum value, then in the 
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rest frame of X2 both xi an d the £ + £ pair are forced to be at rest. The momentum of X2 m the 
laboratory frame is then determined 

i^o = (l + M^/M e+£ -)P i+t - . 

where Pe+e- is the momentum of the dilepton system. The x® can then be combined with b jets to 
reconstruct the decay chain. A clear correlation between the masses of the 6x2 anc ^ bbyQ systems 
was observed allowing both the gluino and sbottom masses to be determined if the mass of Xi was 
assumed. The inferred mass difference rrig — was found to be insensitive to assumed Xi mass. 

In the case of interest here the situation is more complicated. First, there is an extra step in the 
decay chain i.e. g —* bb — ► bbx^b —* bbrf — > bbr + T~Xi- So that even if the events could be selected 
such that the rr invariant mass was at the kinematic limit, Xi would not be at rest in the x.2 res ^ 
frame, and the inferred X2 momenta would not be correct. This was the case at "Point 5" [ffl] where 
the method was applied to the decay chain q — ► qx% — * — ^ Q^fJ- Xi an d nevertheless, a mass 
peak was reconstructed in that case. Second, the momentum of the rr system cannot be measured 
owing to the lost energy from neutrinos. Despite these problems the method is still effective as is 
now demonstrated. We select events with reconstructed rr mass in the range 

40 GeV < m TT < 60 GeV 

and infer the momentum of X2 from the measured momentum P T + T - of the rr system assuming 
the nominal value of M^o . 

P^o = (l + M^o/M T+r -) P T+T - . 

This momentum is then combined with that of two measured b-jets each required to have Pt > 
25 GeV and the mass of the X2P ano ^ X^bb systems computed. Figure || shows the correlation 
rn{x2,b) vs. (m{x^bb) — m{x^2,b)) in a lego plot. The subtracted distribution corresponding to r + r~ — 
r _ r _ — r + r + is used to reduce the background. There is a clear peak in this plot. The projection 
of this plot onto the m{y^bb) — m^ft) axis is shown in Figure which shows a peak at 120 
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Figure 7: Projection of Figure [6| onto the m^x^bb) — m(x2&) axis. 
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Projection of Figure |6| onto the m{x^bb) axis with the requirement that 100 GeV < 
m(x%b) < 140 GeV. 



GeV, somewhat below the true mass difference of 150 GeV. If a selection of events with 120 GeV 
< ?7i(x2&&) — m(x2&) < 140 GeV is made and Figure ^ projected onto the m(%2^) ax ^ s the result 
is shown in Figure ^. A fairly sharp peak results at a value somewhat below the gluino mass of 
540 GeV. This displacement to lower values is due to two effects; jet energy is lost out of the 
clustering cone and carried off by neutrinos is semileptonic bottom and charm decays. We have 
not recalibrated the b— jet energy scale to take account of these effects. 
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Figure 9: Reconstructed rr + jet mass distribution at for light quark jets. 

5 Light Squarks 

We now attempt to find evidence for the decay chain qi — > qx® ~^ Qtt — ► qTTXi- The rates 
are not large due to the small branching ratio for the first step, and we can expect considerable 
combinatorial background from QCD radiation of light quark and gluon jets. The event sample of 
section |3| is used. In addition we require the presence of a non 6-jet with Pt > 25 GeV. Events are 
selected that have reconstructed tau pairs with invarient mass within ±10 GeV of peak in Figure 
and the invariant mass of the tau pair and the jet is formed. This mass distribution is shown in 
Figure |9[ The sign subtracted distribution corresponding to t + t~ — t~t _ — t + t + is used as it 
reduces combinatorical background. There should should be an edge at ~ mq, — ~ 400 GeV. 
The edge is not sharp — two ^ r 's and the Xi are an lost • I n addtion the distribtution is contaminated 
by decays from X3 an d Xi- While this distribution is clealy distinct from that shown above where 
b-jets were used, more work is needed to establish that this could be used to infer information on 
the light squark mass. 



6 Extraction of q~R 

This analysis is based on the fact that <jr — > qx± is dominant, so q~Rq~R pair production gives a 
pair of hard jets and large missing energy. There is no kinematic endpoint, but the Pt of the jets 
provides a measure of the squark mass 0. The following cuts were made: 

• $ T > 200 GeV 

• 2 jets with P T > 150 GeV 

• No other jet with px > 25 GeV 

• Transverse sphericity St > 0.2 
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Figure 10: Transverse momentum distribution for jets passing the selection described in section 
The Standard Model background is shown as the filled histogram. 



• $ T > 0.2M efT 

• No leptons, No b-jets, No tau jets 

The transverse momentum distribution of the leading jets is shown in Figure 0. The error on 
the mass is limited by the systematics of understanding the production dynamics and the SUSY 
backgrounds. Studies of other cases H have shown that this distribution should enable a precision 
of ±50 GeV to be reached; it might be possible to achieve ±25 GeV in a high statistics study. 



7 Dilepton Final states 

While the light gauginos decay almost entirely into r's, the heavy ones can decay via X3,4 ~~ > 
il r F^ — ► Xi 2^ + ^ > giving opposite-sign, same-flavor leptons. The largest combined branching is for 
X% —> ^ ~^ Xi^ + £~ j which gives a dilepton endpoint at 



MI 



\ 



163.2 GeV 



There is of course a large background from leptonic r decays, but this can be cancelled statistically 
by measuring the flavor-subtracted combination e + e~ + /x + pT — e^/U^ as we now demonstrate. 
Events were selected to have two leptons with Pj- > 10 GeV and \rj\ < 2.5 in addition to the jet 



and fx cuts described earlier (see Section y: no tau requirement is applied here). Figure 11 shows 
the distribution in the final state. A clear peak from Z decay is visible that results from 

X3 and X4 decays. The flavor-subtracted combination e + e~ + p + n~ — e^p^ is shown in Figure 12 



and shows an excess extending to ~ 160 GeV. Unlike the distributions involvin tau final states, 
this one can be extrapolated to high luminosity operation which will surely be needed to extract a 
quantitative result from it. 
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Figure 11: Reconstructed /x + /U mass distribution.. The filled histogram shows the Standard Model 
background. 
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Figure 12: Reconstructed + e + e 



fj, e + mass distribution. 



8 Determining SUSY parameters and Conclusion 

The presence of the dijet signal of Section |6| implies that that m g > rriq R . Likewise the failure to 
observe a dilepton peak implies that m e - fl > m^o . These results are used together with the assumed 



■mi 



errors on the measured quantities to fit the model parameters. The values of the errors on m g - 
m g , m,q r and the rr edge are shown in Table |2[ We do not use the information from Figures ^ and 
12 as we have not estimated the quantitative information that they could give. Two fits are shown 
since the sign of fj, cannot be determined. This is expected: a change of conventions can replace 
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Table 2: Results of the fit for the model parameters. The assumed errors in GeV on the measured 
quantities are shown for low and high luminosity with two different assumptions about how well 
the q~R mass can be extracted from Figure [l0|. The fitted values of mo mi/2, t an & an d ^4o are given 
for each case for both signs of \i. The theoretical plus experimental error on the light Higgs mass 
is assumed to be 3 GeV. 
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m 


230 ± 37 


232 ± 26 


230 ± 26 


233 ± 26 


mi/2 


200 ± 14 


196 ± 11 


198 ± 7 


201 ±6 


tan (3 


42 ± 7.3 


42 ±7.1 


45 ±6.2 


45 ±6.1 


A 


0±270 


0±270 


-100 ± 210 


-150 ± 200 



sgn/i with sgn(tan/3), and tan/3 = ±oo are equivalent. 

We assume that the Higgs mass is measured via its decay to two photons. The error on the Higgs 
mass is likely to be dominated by the theoretical uncertainty on the higher order corrections; both 
the one-loop and the dominant two-loop effects have been calculated and are large. The present 
error is probably about ±3 GeV; this might be reduced to ±1 GeV with much more work. The 
ultimate limit comes from the experimental error, about ±0.2 GeV. The effect of reducing this error 
is only apparent in the error of the fitted value of tan (3 whose error is reduced by approximately a 
factor of two if ±1 GeV error on the Higgs mass is used. The table shows various assumptions for 
the errors that might be achieved. The numbers in the first column are conservative and will be 
achieved with the 10 fb _1 of integrated luminosity shown on the figures. The rightmost column is 
an estimate of what might ultimately be achievable. We caution the reader that the measurements 
involving tau's may not be possible at a luminosity of 10 34 cm -2 sec -1 due to pileup effects. 

We can see from the table that, despite the fact that the tau momenta cannot be measured 
directly due to the presence of neutrinos in their decays, we can still expect to infer values of the 
underlying parameters with errors of better than 10%. Of course these errors are considerably 
poorer than those that we expect in cases where taus do not have to be used Q. Our encouraging 
result arises mainly from the very large statistical sample that LHC can produce for the case 
considered. 
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